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We investigate the quantum properties for the codirectional three-mode Kerr non- 
linear coupler. We investigate single-, two- and three-mode quadrature squeezing, 
Wigner function and purity. We prove that this device can provide richer nonclas- 
0^ ■ sical effects than those produced by the conventional coupler, i.e. the two-mode 

| Kerr coupler. We show that it can provide squeezing and the quadrature squeez- 

ing exhibiting leaf-revival-collapse phenomenon in dependence on the values of the 
interaction parameters. In contrast to the conventional Kerr coupler two different 



forms of cat states can be simultaneously generated in the waveguides. We deduce 



conditions required for the complete disentanglement between the components of the 
system. 

PACS numbers: 42.50Dv,42.60.Gd 
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I. INTRODUCTION 



The optical coupler is a device composed of two (or more) waveguides, which areplaced 
close enough to allow exchanging energy between waveguides via evanescent waves R- e_ 



2 



cently, this device has attracted much attention for several reasons. The progress in the 
optics communication and quantum computing networks requires data transmission 2||. 

n n 

This simple device has potential applications in all-optical switching [3|, |4J. Furthermore, 
it provides electromagnetic fields with an exceptionally wide range of nonclassical effects. 
Most importantly this device has been implemented p, |6| and applied in many experimen- 
tal approaches, e.g. in picosecond switching induced by saturable absorption [7], optical 
multi-mode interference devices based on self-imaging [8J and photonic bandgap structures 
in planar nonlinear waveguides 9]. Also the generation of correlated photons in controlled 
spatial modes by downconversion in nonlinear waveguides has been discussed in flol |. 

Quantum mechanically, there are different types of directional couplers. For instance, 
symmetric coupler lll| (linear (or nonlinear) processes are involved in both the waveguides), 
asymmetric coupler [l2j (at least one of the waveguides possesses different nonlinearity than 



13, 1181, 
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20, 



review papers 
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23) 



the others), Raman-Brillouin coupler jl3], bandgap coupler |l4j and Kerr coupler [15| , 
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is 



221 ]. For more details about their properties the reader can consult the 



Kerr nonlinear coupler (KNC) has taken a considerable interest in the literature since the 
third-order nonlinearity provides an effective mechanism for generating nonclassical effects in 
the electromagnetic waves caused by the processes of self-phase and cross-phase modulations. 
Also this is related to potential advantage of the possible observation of the large values of 
third-order optical nonlinearities in the organic polymers jfj . In this regard the generation of 
continous variable Einstein- Podolsky- Rosen entanglement via Kerr nonlinearity in an optical 
fiber [24j as well as the generation of spatial soliton arrays in a planar Kerr waveguide from 
seeded spontaneous parametric down conversion [25j are achieved. Several papers have 
been devoted to the quantum properties of the KNC in the framework of rotating-wave 
approximation by neglecting the rotational terms in such a way that a closed form solution 



is obtained 
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[19I . [2^]. In this case KNC has provided many of interesting 



effects such as revival-collapse phenomenon (RCP) in the mean-photon numbers, squeezing 
of vacuum fluctuations, sub-Poissonian statistics in single as well as in the compound modes 

n 

|15l ] . The phase distribution of KNC has been investigated showing that the phase-difference 
evolution is closely connected with the energy exchange between waveguides and the RCP 
in the mean-photon numbers is due to the bifurcation of the phase-difference probability 
distribution [lTj] . Furthermore, the geometry of the waveguides has been considered via 
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varying linear coupling coefficients in the codirectional KNC 16| and contradirectional KNC 
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20| | . In these cases it has been shown that there is a possibility to control the switching 



characteristics and principal squeezing effect by adjusting the shape of the waveguides. Quite 
recently, we have investigated the single-mode quantum properties of the codirectional Kerr 
nonlinear coupler when the frequency mismatch is involved and a condition for obtaining 
an exact solution for the equations of motion is fulfilled 
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22 . For this case we have 



shown that the mean-photon numbers exhibit oscillatory behaviour rather than revivals and 
collapses. Additionally we have proved that the Schrodinger-cat states, in particular, Yurke- 
Stoler cat states (YSCS) [26( can be generated. Also the higher-order squeezing has been 
investigated [22]. 

Till now the KNC has been treated as a two-mode device. In the present paper we give 
for the first time-as far as we know-the three-mode version of this device. The motivation 
of developing such device is that the three-mode KNC can provide nonclassical effects richer 
than those obtained from the conventional coupler (i.e. the two-mode version), as we shall 
show throughout the paper. As is well known that the basic efforts in the quantum optics 
is to enhance the nonclassical effects. Additionally, we show that the quadrature squeezing 
exhibits leaf-revival-collapse phenomenon in dependence on the values of the interaction 
parameters. In this phenomenon the revival patterns provide leaf shapes and between two 
revival patterns short collapse period occurs. It is worth reminding that for the standard 
revival-collapse phenomenon the revival patterns have ellipsoid shape. We proceed that, in 
contrast to the conventional Kerr coupler two different forms of cat states can be simultane- 
ously generated in the waveguides. The investigation of the three-mode KNC will be given 
in the following order: In section 2 we give the Hamiltonian for the system and the solution 
for the equations of motion. Also we derive the expectation values for different moments 
of operators, which will be used in the paper. In section 3 we investigate the quadrature 
squeezing. In section 4 we discuss the evolution of the Wigner function and the purity for 
the single-mode case. In section 5 we give the main conclusions from the results. 



II. MODEL FORMALISM AND DYNAMICAL SOLUTION 

In this section we give the Hamiltonian for the system and derive the solutions for its 
equations of motion. 
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The Hamiltonian controlling the three-mode codirectional Kerr nonlinear coupler can be 
represented as 

3 



H 

h 

3=1 

(1) 



+Ai[dia2 exp(zAit) + djd2exp(— iAit)} + A2[did3 exp(zA 2 t) + a[as exp(— zA 2 t)], 

where A± = u)i — lu 2 , and A 2 = ui\ — c^. The waves are designated by di (fundamental 
and/or first), a>2 (second) and 03 (third) modes with frequencies ct>i,o>2 and respectively. 
The coupling constants Xj an d Xj are proportional to the third-order susceptibility x 
and are responsible correspondingly for the self-action and cross-action processes of the jth 
mode. The linear coupling between the waveguides are represented by Ai and A2, which we 
assume to be real. Hamiltonian (CQ) gives a generalization to several models discussed in the 
literatures earlier by controlling the values of the interaction parameters such as anhramonic 
oscillator 



271 ]. up-conversion process [28( and two-mode KNC |2l|. The scheme describing 



Hamiltonian ([I]) is shown in Fig. 1 and can be explained as follows. Two waveguides are 
operating by Kerr-like nonlinear processes. In the first waveguide the fundamental mode 
di propagates, while in the second waveguide the second d 2 and third d 3 modes propagate. 
The interaction between the fundamental mode and the second-third modes occurs via the 
evanescent waves. We have assumed that the all waves are propagating with the same 
velocity v , hence the time t and travelled distance z are related by z = vt. Outgoing fields 
can be detected as single or compound modes by means of homodyne, photocounting or 
coincidence detection in the standard way. Throughout the investigation of the system we 
do not consider dissipation, which generally decreases the amount of nonclassical effects. 

We have to comment that the linear coupling coefficient is from the wave equation equal 
to X/v = co 2 x^ /2kc 2 , whereas the nonlinear coupling coefficients are x/ v = uj 2 x^ /%kc 2 (v 
being the speed of light in the medium, u being frequency of the travelled wave and k is its 
wave number). Taking into account that values of quadratic susceptibilities lie in an interval 
10~ 14 — 10~ w m/V and values of cubic susceptibilities in an interval 10~ 17 — 10~ 13 m 2 /V 2 for 
various nonlinear materials, then in the optical region the values of X/v lie from 10 4 to 
10 7 m _1 and values of x/ v from 10 to 10 5 m -1 . So the effects predicted should be expected 
to be observable in samples of the length of centimeters and in nanosecond scales. 

Now introducing new operators Aj = aj exp(iu)jt), j = 1, 2, 3 we can write the Heisenberg 
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FIG. 1: Scheme of realization of interaction (1). 



equations for ([I]) as 



^ = -i(2xiilii + Xii|i 2 + x 2 i^i 3 )ii - iAii 2 - zA 2 i 3 , 



^ = -i(2 X2 AlA 2 + xA^i + X 3 AlA 3 )A 2 - i\ 1 A 1 , 



dA 3 
dt 



i(2 X3 AlA 3 + x 2 ^i^i + x 3 A\A2)A 3 - i\ 2 A 1 . 



(2) 



To solve ([2]) exactly we assume the processes of the self-action and cross-action compensate 
each other in the evolution of the system. This can be expressed as X — Xi — X2 = Xzi X — 
Xi = X2 = X3 an d 2% = X- m i^-ct, the cross-spectral coupling between various waveguides 
may be comparable with the self-coupling provided that the surfaces of the waveguides are of 
high quality (roughness is much less than the wavelength) to avoid a substantional reduction 
of amplitudes by evanescent damping. As follows from the values for quadratic and cubic 
susceptibilites both the linear and nonlinear couplings may be simultaneously significant for 
sufficiently strong fields. Under these conditions one can easily prove that 



A\A ± + A\A 2 + AlA 



3^3 



N 



(3) 



is a constant of motion. Based on this fact the system of equations (T2J) can be modified via 
the substitution Bj = exp(2ixNt)Aj, j = 1, 2, 3 and takes the form 
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at 



3- 



(4) 

The system (j3J) can be easily solved, e.g., using Laplace transformation, and the general 
solution is 

Ai(t) = exp(-2zxiV£)jai(0) cos fit - zJ(Aia 2 (0) + A 2 a 3 (0)) sin /it}, 

i 2 (t) = exp(-2i X iVt){[l - 2$ sin 2 (f )]a 2 (0) - sin 2 (f ) - sin(/it)}, (5) 

i 3 (t) = exp(-2* X iVi){[l - 2p sin 2 (f )]a 3 (0) - 2^M sin 2 (f ) - sin(^)}, 



where /i = a/A| + A§. There are several facts which can be extracted from the structure of 
the Hamiltonian ([1]), the solution ([5]) and Fig. 1. When a 3 (0) <-> a 2 (0) and Ai <-> A 2 one 
obtains A 2 (t) <-> A 3 (t). The nature of the coupler, i.e. the switching of energy between 
waveguides, manifests itself by periodic functions in ([5]), whereas the Kerr nonlinearities in 
the waveguides are described by the nonlinear (quadratic) phase, which plays an essential 
role in generating the nonclassical effects. On the other hand, the fundamental mode a\ 
can provide richer nonclassical effects than those produced by the second and third modes 



as well as the conventional coupler [2l|, |22j (see Fig. 1). This is resulting from the coupler 
mechanism, which switches the energy in the second-third mode (second waveguide) jointly 
to the fundamental (first waveguide) and vice versa. 

On the other hand, the nature of the field quantization is evident where one can easily 
prove from ([5]) that 

i,(t),4(t)] =6 jd ,, (6) 

where 6jj> is the Kronecker delta and j,j' = 1, 2, 3. In the derivation of ([6]) one has to use 
the identities 



exp(-2ixA r t)a i (0)exp(2zxA r i) = a j (Q)exp(2ixt), 

(7) 

exp(-2ixNt)a](0) exp(2ixNt) = o}(0) exp(-2z'x£). 
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Furthermore, the mean-photon number for the fundamental and second modes are: 
A[{t)Ai{t) = ai(0)al(0) cos 2 (fit) 



+ 



sin (pt) 



A?oS(0)o 2 (0) + A 2 4(0)a 3 (0) + A X A 2 (a 2 (0)a 3 (0) + a 2 (0)a 3 (0) 
+i[A 1 (a 1 (0)aJ(0) - al(0)a 2 (0)) + A 2 (a 1 (0)a|(0) - a\(0)a 3 (0))] a ^ 

A\(t)A 2 {t) = 4^(0)^(0) sinV) + [1 - 2$aJ(0)a 2 (0) sin 2 (f )] 2 + 4^4(0)a 3 (0) sin 4 (f 



-2^[1 - 2$ sin 2 (f )] sin 2 (f ) [a 2 (0)a 3 (0) + a 2 (0)4(0)] - ^[1 - 2$ sin 2 (f )] sin(^) 

x [ax (0)4(0) - 4(0)a 2 (0)] + ^sin 2 (f ) sin(//t) ^(0)4(0) - 4(0)a 3 (0)]. 
From ([S]) the correlation between modes and switching mechanism, i.e. when one photon is 
created in the one of the modes the other is annihilated in the other mode, are quite obvious 

3 „ m 

We close this section by calculating the moments (Y[ A™ 3 A™ 3 ) when the modes are 
initially prepared in the coherent states \a.\, a 2 , a 3 ) and aj are real. From §5§ one can easily 
deduce that 

3 

<n^ tna r> = ar(*)^ ii (*)^ na (*)s? 9 («)fir"(<)ss"(*)« fci *M<z h2 - 1)], (9) 

3=1 

where nj,rrij are positive integers, 

e = |«i| 2 + |a 2 | 2 + |a 3 | 2 , z = exp(-2ix^), 

3 

hi — I ^2[mj(mj - 1) - nj(rij - 1)] + m l {m 2 + m 3 ) + m 2 m 3 - n x (n 2 + n 3 ) - n 2 n 3 , (10) 



/i 2 = mi + m 2 + m 3 - n x - n 2 - n 3 
and 

a 1 (t) = a x cos(/it) - z[Aio; 2 + A 2 a 3 ]^^ = a lx (t) + ia ly (t), 

a 2 (t) = a 2 [l - ^sin 2 (i^)] - sin 2 (f ) - z^f sin(^) = ^ (t ) + ia 2y (t ) , (11) 

a 3 (t) = a 3 [l - ^sin 2 ^)] - sin 2 (f ) - i*f> sin(^) = a^t) + ia 3y (t). 
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Expression (Q will be frequently used in the paper. It is obvious that when \t = inn 
and m is integer, i.e. z — 1, this expression reduces to that for the coherent light with 
amplitudes ay (t) indicating that the system produces coherent light periodically. Also from 
this expression one can easily check that (A*™' 1 A™ 1 } = (A^Aj)™* . This means that the 
system is Poissonian, i.e. it cannot exhibit sub-Poissonian statistics, and the mean-photon 
numbers cannot provide RCP. Expression (Q reduces to that of the two- mode KNC 21 ] 
when A 2 = and a 3 = 0. Now we use the calculations given here to investigate quadrature 
squeezing, Wigner function and purity in the following sections. 

III. QUADRATURE SQUEEZING 

Squeezing is a pure nonclassical phenomenon without classical analog. Squeezed light 
has less noise in one of the field quadrature than the vacuum level and an excess of noise 
in the other quadrature such that the uncertainty principle is satisfied. Squeezed light can 
be measured in the homo dyne detection. This light has a lot of applications, e.g., in optical 



communication networks [29(, quantum information [30|] and high precision measurements. 
In this section we investigate the single- and two-mode squeezing and shed the light briefly 
on the evolution of the three-mode squeezing. In doing so we define two quadratures X n 
and Y n as 

1 n 1 n 

Xn = ~ Y^lMt) + Yn = ~ Y^Mt) ~ (12) 

3=1 1 3=1 

where n takes on values 1, 2 and 3 associated with the single-, two- and three-mode squeezing, 
respectively. Quadratures (fl2"j) satisfy the following commutation rule: 

[XnX} = ^ (13) 

where C n is a c-number and takes on the value C\ — 1, C2 = 2, C3 = 3. The uncertainty 
relation associated with the commutation rule ( fl3|) is 

<(AX n ) 2 )((AK„) 2 >>^, (14) 

where ((AX n ) 2 ) = (A 2 ) — (A n ) 2 and similar form can be given for ((AY n ) 2 ). The system 
is said to be squeezed in the A n -quadrature if 

S n = 4((AA n (t)) 2 ) - \C n \ < 0. (15) 
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The equality sign in (IT51) holds for minimum-uncertainty states. Similar definition can be 
given for the F n -quadrature (defining a Q n -factor). Based on Q and ffl~5l) we investigate 
the single- and two-mode squeezing in details and comment on the three-mode squeezing. 
As we deal with the Kerr media we can expect squeezing in the guided modes in all these 
types. 

A. Single-mode squeezing 

The squeezing factors for the jth mode can be straightforwardly evaluated as 
S[ j \t) = 2<xj x (t) + 2<xj y (t) + 2^[a] x (t) - a 2 y {t)} cosOi + 2a jx (t)a jy (t) sin6i} exp[-2esin 2 (2xt)], 

—4 [ocj X (t) cos6 2 + otjyit) sin 2 ] 2 exp[— 4e sin 2 (xt)], 

(16) 

Q[ j \t) = 2a 2 x {t) + 2a) y (t) - 2^[a 2 jx (t) - a 2 y {t)} cosGi + 2a jx {t)a jy {t) sinGi} exp[-2e sin 2 (2xt)] 

—4 [otj X {t) sin 2 — ocj y {t) cos 6 2 ] 2 exp[— 4e sin 2 (xt)], 
where 

0! = 2 X t + esin(4xt), 6 2 = esm{2 X t). (17) 

The subscript 1 means single-mode squeezing and the superscript j denotes the mode under 
consideration. Now we give some analytical results showing that the system can provide 
single-mode squeezing depending on the values of the interaction parameters. For instance, 
when («i, a 2 , 03) = (a, 0, 0) and xt — fnir/2, m is odd integer, expressions ( ITBl) give for the 
fundamental and second modes the following: 

S { ^\t) = -4a 2 cos 2 ( / ut)exp(-4e), Q^(t) = 4a 2 cos 2 (/it), 

(18) 

S[ 2 \t) = 4^f sinV), Qf\t) = -A$f sin 2 (/it) exp(-4e). 

It is worth reminding that the expression associated with the third mode can be obtained 
from that of the second mode via the transformation Ai <-> A 2 and a 2 <-> 0:3. From ( fl8l) 
the complementarity between the squeezing factors in the two waveguides is obvious, i.e. 
maximum squeezing in the first waveguide is accompanied by minimum squeezing in the 
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FIG. 2: Evolution of the single-mode squeezing factors for x = 0.5s -1 , (02, 03) = (0.3, 0.3) and for 
(a) (ai,Ai,A2) = (0, Is -1 , Is -1 ) (solid curve for Q^'(t), short-dashed curve for s[ 2 \t) and long- 
dashed curve is the squeezing factor Q^\t) for the two-mode KNC when (a\, 02, 03, Ai, A2) = 
(0,0.3V2s- 1 ,0,ls- 1 ,0s- 1 )), and (ai,Ai,A 2 ) = (0.3, Is" 1 , 50s" 1 ) for (b) s[ 1] (t) and (c) s[ 2) {t). 

second waveguide and vice versa. Maximum squeezing occurs in the fundamental mode 
when /it = mV and ml is integer. In this case, (t) can be expressed as 

Si (t = mV//i) = — rj exp(— rj), r] = 4a 2 . (19) 

By evaluating the extreme values for ( fT9l) one finds that the maximum squeezing in S± (t) 
is 37% when a = 0.5. Also form (TIBl) it can be proved that the maximum squeezing 
involved in S[ (t) is the half of that of S^^i), i.e. 18.5%. It is worth reminding that for 
the single-mode case C\ = 1, i.e. S±\t) is normalized. Now we clarify the role of the 
switching mechanism in the behaviour of the coupler since the Kerr-like media can generate 
nonclassical effects regardless of entanglement between modes. This can be understood by 
assuming that Ai = A2 = 0, i.e. the linear interaction between waveguides is neglected, and 
c*i = 0. In this case the fundamental mode cannot generate nonclassical effects, however, 
in the framework of coupler the fundamental mode exhibits nonclassical effects. In other 
words, entanglement between different components of the system can generate as well as 
amplify the nonclassical effects (if they exist). Information about this situation is shown 
in Fig. 2(a) for the given values of the parameters. It is worth mentioning that for Qi (t) 
of the two-mode KNC (as well as S2 (t) in Fig. 3(a)) the initial total photon number has 
been considered equal to that of the three-mode KNC for significant comparison. In this 
case we have noted that the squeezing occurs only in the one of the quadratures, i.e. in 
Y n (X n )-quadrature for the fundamental (second) modes. Also squeezing generated in the 
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fundamental mode is much greater than that produced in the second mode as well as in the 
single-mode of the two-mode KNC Influence of the values of the linear coupling constants 
Xj, i.e. the intensities of energy switching between the waveguides, on the evolution of the 
quadrature squeezing is shown in Figs. 2(b) and (c) for Sf^(t) and S[ (t), respectively. It 
is worth mentioning that Qi (t) provides a behaviour typical to that shown in Fig. 2(b). 
From Fig. 2(b) one can see that additionally to the occurrence of the nonclassical squeezing 
the s[ (t) exhibits leaf-revival-collapse phenomenon. Actually, there are two requirements 
for occurring such phenomenon, which are: (i) The linear and nonlinear couplings in the 
Hamiltonian (CQ) should contribute significantly to the dynamics of the system, (ii) The 
energy exchange between the mode under consideration with at least one of the modes in the 
other waveguide is very strong. With this in mind we can understand why S[ 2 ^ (t) provides 
squeezing rather than RCP. This can be realized easily when A 2 >> Ai, i.e. Ai/A 2 — 0, 
a 2 (t) = a and hence the time evolution of S^\t) basically depends on the contribution of 
the Kerr-like medium, which alone cannot produce RCP. Generally, the behaviour of the 
system for this case can be explained as follows. As A 2 >> Ai, i.e. the energy switching 
mainly occurs between the fundamental and third modes, which behave as a two-mode KNC 
and the second mode practically does not affect the behaviour. Consequently, the second 
mode behaves as an independent mode evolving in a Kerr medium. In conclusion, when Ax 
and A 2 are large or Ai >> A 2 the squeezing factors for the fundamental and second modes 
provide RCP. On the other hand, the occurrence of the RCP in the squeezing factors of 
the fundamental mode can be explained in the following sense 2jJ. The squeezing factors 
(fT6l) include two forms of periodic functions: one is coming from the self-cross-nonlinear 
interaction part of (pQ), in particular, the exponential function whose period is n/x> an d 
the other is arising from the linear- interaction part whose period is ir/fi. When the linear 
interaction between waveguides is very strong-either between the fundamental- second or 
fundamental-third or both-the period of the energy exchange between waveguides decreases, 
i.e. many oscillations occur, till the interaction time becomes t = tt/x, at this moment the 
field is trapped instantaneously by nonlinearity in the waveguides and the squeezing factors 
show collapse. The oscillations of the RCP are caused by the linear coupling Xj and the 
envelope of the revivals is caused by the nonlinearity in the system. As the interaction 
proceeds the phenomenon is periodically repeated. 

We have to stress that the evolution of the squeezing factors is sensitive to the values of 
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the intensities, however, in the strong-intensity regime it can provide nonclassical squeezing 
instantaneously greater than that present in Figs. 2, but the behaviour becomes rather 
complicated. We draw the attention to this point in the following section. 

B. Two-mode squeezing 

In this part we study the two-mode squeezing in which the correlation between modes 
starts to play a role. Also we compare the behaviour here with the corresponding one for the 
two-mode KNC We restrict the analysis to the fundamental-second (between waveguides) 
and second-third modes (in the same waveguide). The fundamental-second mode squeezing 
factors can be straightforwardly evaluated as 

S^ 2 \t) = I[5i 1} (t) + S?\t)) + 2{{a lx (t)a 2x (t) - a ly (t)a 2y {t)) cos9 2 

+[a 2x (t)ai y (t) + a lx (t)a 2y (t)} sin6 2 | exp[-2e sin 2 (2xt)] + 2[a lx (t)a 2x (t) + a ly (t)a 2y (t)} 

—4 [ai x (t) cos 0i + aci y {t) sin X ] [a 2x (t) cos ©i + oc 2y (t) sin X ] exp[— 4e sin 2 (xt)], 

(20) 

Q { 2' 2) (t) = l[Q { i\t) + Q?\t)] - 2{[a lx (t)a 2x (t) - a ly (t)a 2y (t)} cos0 2 

+ [a 2x (t)ai y (t) + a lx (t)a 2y (t)] sin6 2 | exp[-2e sin 2 (2xt)] + 2[a lx (t)a 2x (t) + a ly (t)a 2y (t)) 

—4 \o>i y {t) cosGi — ai x (t) sin0i] [a 2y {t) cosOi — a 2x {t) sin 0i] exp[— 4e sin 2 (x^)], 

where S[ j) and Q ( f are given by (|T6l) and Bj have the same meaning as given above. The 
subscript 2 stands for two-mode squeezing and the superscript (1,2) denotes fundamental- 
second squeezing factor. The second-third squeezing factors can be easily obtain from (|20|) 
by simply replacing 1 by 3. As we did above we analyse the case (ai, a 2 , a 3 ) = (a, 0, 0) and 
\t = mir/2, m is odd integer, for which expressions ( |20l) reduce to 

4 1 ' 2) (t) = 2a 2 [^|sin 2 (^)-cos(/it)exp(-4 e )], 
Q { z' 2 \t) = 2a 2 [cos(/it) - ^|sin 2 (/it)exp(-4e)]. 

Squeezing can be periodically realized in the two quadratures for particular values of in- 
teraction time. The two-mode squeezing between the second and third modes, i.e. modes 
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° 25 0.0 1.0 2.0 3.0 4.0 5.0 6.0 0.0 2.0 4.0 6.0 8.0 10.0 12.0 

t t 

FIG. 3: Evolution of the two-mode squeezing factors for x = 0.5s -1 , («i, a>2, 03) = (0,0.3,0.3) 
and for (a) Ai = Ai = Is -1 (long-dashed curve for sj 1 ' (t), solid curve for S^' 3 \t), the 
short-dashed curve is given for S^i) factor of the two-mode KNC when (a±, 02, «3, Ai, A2) = 
(0,0.3^/2,0, Is" 1 , 0s -1 )) and the straight line is given for showing the squeezing bound), (b) 
S^ 3 \t) for Ai = 50s" 1 , A 2 = Is" 1 . 

evolve in the same waveguide, for this special case is 



4 2 ' 3) (t)=2a 2 (^sinV), 

(22) 

Q?' 3) (£) = -2a 2 ^^sin 2 (/it)exp(-4 e ). 

Squeezing can occur in the second quadrature only similar to the single-mode case. Now 
we draw the attention to the behaviour of the general forms ( 1201 . which are presented in 
Figs. 3(a) and (b) for given values of the interaction parameters. From Fig. 3(a) one can 

(1 2") 

observe that S 2 ' (t) and S 2 (t) exhibit nonclassical squeezing periodically. Complementary 

(\ 2) 

behaviours have been noticed for Q 2 ' (t) and Q2(t). Nevertheless, for modes evolving in 

the same waveguide where switching mechanism does not exist directly, squeezing occurs-for 

(2 3) 

these values of interaction parameters-only in one of the two quadratures, i.e. S 2 ' (t), as 
indicated by the solid curve. We have to stress that the three-mode KNC provides amount 
of squeezing in the framework of two-mode squeezing greater than that of the two-mode 
KNC (compare the solid curve with the short-dashed curve). The action of the coupling 
constant is shown in Fig. 3(b) for sif' 3 \t). From this figure it is obvious that squeezing 
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(or minimum-uncertainty light) always occurs and the amount of squeezing is greater than 
that obtained for weak-coupling case (compare with the solid curve in Fig. 3(a)). This 
means that the linear mechanism in the system can play a significant role in amplifying the 
nonclassical effects. RCP is remarkable, which always occurs provided that Ai or A2 or both 
are large. Also for these cases we have noted that the two-mode squeezing factors between 
the waveguides exhibit RCP somewhat similar to that shown in Fig. 2(b). 

Finally, for the three-mode squeezing factors we have noted that squeezing periodically 
occurs in the two quadratures similar to that of the two-mode case, in particular, the 
fundamental-second mode factors. Moreover, the amount of squeezing achieved in the three- 
mode factors are greater than those in the single-mode and two-mode factors. Also intense 
switching between the waveguides provide RCP in the evolution of the three-mode squeezing 
factors. 



IV. QUASIPROB ABILITY DISTRIBUTION FUNCTION 

Quasiprobabihty distribution functions, namely, Wigner W-, Husimi Q- and Glauber 
P-functions |3jj, provide all physical information about the quantum mechanical systems. 
Investigation of these functions for the quantum mechanical systems is one of the main 
topics in quantum optics. Actually, among all quasiprobabihty functions W function takes 
a considerable interest since it has been experimentally realized by different techniques, e.g. 
in homodyne tomography [32], photon counting experiment (3] and in experiments with 



trapped ions 3J], and it is sensitive to the interference in phase space. In this section we 
investigate the single-mode quasiprobabihty distribution functions, in particular, W function 
for the system under consideration. As the derivations of these functions are similar with 

n 

those in [21|, we write down only the explicit formulae for them. Just here we give only 
outline of the derivation. The solution (jSJ) includes exponential quadratic operator phase, 
which causes difficulties in evaluating the qu asip robabihty functions in the standard way. 
Thus one has to use the technique given in [35(, which basically depends on the formula 
(Q, to achieve the goal. It is worth reminding that the modes are initially prepared in the 
coherent state |ai, c*2, CK3). Furthermore, we extend the investigation to include the purity 
for the single-mode case, which can be evaluated by means of the symmetric characteristic 



function. Now the jth mode symmetric characteristic function is 21] 
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FIG. 4: The W function for the single-mode case for \ = 0.5s ,i = ns and for ay = 0.9, Ai = 
A2 = l/(2\/2s _1 ) (a) (fundamental mode), (b) (second mode); and a\ = 2, 02 = «3 = and 
Ai = A2 = \j\f2s~ 1 (c) (second mode). 



^(C,t) = exp(^|C| 2 ) £ 



ni,n2=0 



xz [^(n 2 -l)-^(m-l)] exp [ e ( z n2-ni _ 

Also the PF function for the jth mode takes the form 211 ]: 



(23) 



^■(/M) = ~exp(-2 



ni,n.2=0 



(24) 



xz 1 2 



[^(na-l)-^(m-l) 



exp[e(z" 2 - ni - l)]L"2-"i(2|/?| 2 ), 

where L™(.) is the associated Laguerre polynomial of order n and (3 = x + iy. One can 
easily verify that when yt — m' 71 and m' is integer, the system produces coherent light with 
amplitude dj(t) 21[. Furthermore, for t\ = [m! + l/2)7r and m' is integer, expression ([24 
gives that for the cat states (the derivation is given in the appendix of |2l|) 



(25) 



Wj((3, t) = ijexp (-2|/5 - iaj(t)\ 2 ) + exp (-2|/3 + ia 3 (t)\ 2 ) 
+2exp ]-2(|/3| 2 + Dj-)] sin (2((3a*(t) + /3*a,-(t)))}, 



where 



ZA, = e 



l«,(t)| 2 . 



(26) 



Roughly speaking, this form tends to that of the YSCS 26] in the jth mode when Dj ~ 
and dj(t) 7^ 0. In contrast to the two- mode KNC, the fundamental and second (or third) 



16 



modes of the three-mode KNC can provide different types of cat states for certain choice 
of the interaction parameters. Now we give some analytical facts related to Dj. There are 
three different special cases, which can be discussed for ( 1251) based on the values of ctj where 
we assume Ai = A2. 

(i) When <x,- = a, j = 1, 2, 3 for the fundamental and second modes, we obtain 

D 1 = a 2 [l + cos 2 (/it)], 

(27) 

D 2 = a 2 [2 + ± sin 2 (/it)]. 

Expressions f l2T|) show that D 2 > D\ and Dj ~ only when a < 1. This means that the cat 
states can be generated in the first waveguide faster than that in the second waveguide. As 



a is relatively small the microscopic cat states [35j] can be generated in the two waveguides 
simultaneously. Trivially, statistical-mixture coherent states can be generated for this case 
in dependence on the values of the interaction parameters, 
(ii) For («i, 0:2,03) = (a, 0,0) we obtain 



D\ = a 2 sin 2 (/it), 

(28) 

D 2 = a 2 [l - isin 2 ^)]. 

Expressions (T28]) indicate that when the fundamental mode provides YSCS in a macroscop- 
ical form a > 1, say, when /it = tt, the second (or the third) produces statistical- mixture 
coherent states, however, when a < 1, YSCS in the microscopic forms are simultaneously 
generated in the two waveguides. This shows the crucial role for the intensities of the initial 
light. 

(iii) For (01,02,^3) = (0, a,0) we obtain 



D 1 = a 2 [l - |sin 2 (/it)], 

(29) 

J D 2 = 2o 2 [l-|sin 2 (ft)]sin 2 (ft). 

This case represents the inverse of the case (ii). The above cases lead to an important fact: 
In order to obtain macroscopically distinguishable YSCS from a particular mode it should 
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be initially prepared in coherent state with large intensity. Information on the above cases 
is shown in Figs. 4 for the fundamental and second modes for given values of interaction 
parameters. The values of the linear constants Xj have been chosen to minimize D\. In 
Fig. 4(a) for the fundamental mode one can observe the two Gaussian peaks and inverted 
negative peak inbetween indicating the occurrence of the interference in phase space. Fig. 
4(b) (for the second mode) provides information on the statistical-mixture coherent states, 
i.e. it shows two-peak forms without interference in phase space. We have to stress that in 
Figs. 4(a) and (b) the states generated are close to those of the microscopic regime. Fig. 
4(c) gives the well-known shape for the W function of the macroscopically distinguishable 
YSCS. In this case the fundamental mode collapses to the state: 



The W function for the second mode corresponding to Fig. 4(c) is that of the vacuum state, 
where 6:2 (t) — 0. Actually, we have noted that such behaviours of the fundamental and 
second modes-related to Fig. 4(c) and W function of vacuum light-can be interchanged 



In Figs. 5 we illustrate the general case, when the above special cases are no longer 
applicable. In these figures we consider t = 12.38s, where squeezing is well pronounced in 
the behaviours of the single-mode case for the fundamental and second modes. In these 
figures the nonclassical effects are presented as asymmetric multi-peak structure indicating 
generation of the multi-component cat states. Also Fig. 5(a) includes additionally nega- 
tive values showing that the fundamental mode gives amount of the nonclassical behaviour 
greater than that for the second (or third) mode. This agrees with the explanation given in 
section 2. On the other hand, when Ai >> A 2 these negative values disappear and the W 
functions of the fundamental and second modes are almost the same (see Fig. 5(c)). It is 
worthwhile mentioning that the corresponding W function for the two-mode KNC gives two 
symmetric peaks without involving any negative values, which is insensitive to the values of 
Ai(> 1). In conclusion, the system can generate different types of YSCS depending on the 
values of the interaction parameters. 

Entanglement characterizes intrinsically quantum-mechanical correlations between quan- 
tum systems. In this regard quantum entanglement is the basic resource required to im- 
plement several of the most important processes studied by quantum information theory. 




(30) 



when one uses t = 7rs, cki = a 3 = 0, a 2 = 2, Ai = A 2 = \/2s 
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FIG. 5: The W function for x = 0.5s" 1 , t = 12.38s, aj = 1,A X = A 2 = ls _i and for (a) the 
fundamental mode, (b) second mode. Figure (c) is given for the fundamental mode for the case 
Ai = 50s -1 , A2 = Is -1 , the values of the others parameters are the same as for (a). 

Thus, we discuss the single-mode degree of purity Tr p?(t) for the system under considera- 
tion, which gives globally information about the entanglement. Purity can be determined by 

the knowledge of quantum state of the system, which in turn can be obtained by quantum 
I I 

tomography [32|. Also it is worth mentioning that in 36] entanglement has been studied 
for the two-mode KNC driven by external classical field showing that the system is able to 



generate Bell-like states. Also in 



371 | the purity and the relative entropy are investigated for 



the pumped dissipative nonlinear oscillator including Kerr nonlinearity. The purity can be 
evaluated by means of the symmetric characteristic function through the relation 



1 



7T 



(31) 



When Trp|(t) = 1 the mode under consideration is in a pure state, i.e. the mode under 
consideration is disentangled from the rest of the system, while when Trp^(t) < 1 the mode is 
in a mixed state and consequently it is entangled with the rest of the system. The subsystems 
are most entangled when their reduced density matrices are maximally mixed. Substituting 
into (13"Tj) and after lengthy but straightforward calculation we arrive at 



Trpj(t) = £ jniimi!] 

ni,mi=0 



(ni+mi)!(-l) m i+"i|a.,(t) | 2 ( n l+ m i> 

(32) 



+ 2 £ £ {J ^SUSS3^S^ ex P[~ 4esin2 ((^2 - m^xt)] cos^(ni,mi,m 2 ), 

mi>m2 ni=0 

where 

■0(ni, 777,1,777,2) = 2xt(m 2 - m 1 )(m 2 - 77,1). (33) 



19 




FIG. 6: The evolution of the purity Tr/5|(i) for \ = 0.5s \ol\ = a 2 = 0.3 and Ai = A2 = Is 1 for 
the fundamental mode (a) and the second mode (b). 

It is evident that fl32l) is dependent on the mean-photon number for the mode under consid- 
eration. Thus for fixed value of x an d xt 7^ 17171 '1 m is integer, the evolution of the purity 
reflects well the switching mechanism of the coupler provided that at least one of the initial 
intensities is large. In relation to Figs. 6, i.e. Ai = A2, a 2 = «3 = ct, the mean-photon 
numbers for the fundamental and second modes take the forms 

a 2 

|ai(t)| 2 = a\ cos 2 (/it) + 2a 2 sin 2 (/it ) , | a 2 (t) | 2 = a 2 cos 2 (/it) + y sin 2 (/it). (34) 

When aj are small, the mean-photon numbers of the fundamental and second modes would 
be relatively small and hence the evolution of the purity will be close to that of the initial 
case (see Figs. 6(a) and (b) for given values of the parameters). Also from Figs. 6 we note 
that the fundamental and second modes tend to partial pure states, i.e. close to pure states, 
when their mean-photon numbers provide their initial values. In this case the components of 
the system are approximately disentangled. Additionally, the fundamental mode only tends 
to almost pure state when the energies in the two waveguides are completely interchanged. 
More illustratively, from these figures one can observe that when t = 0s the two modes are 
in pure states. When t increases and a\ is large, entanglement between modes through the 
interaction ([I]) is achieved and both the fundamental and second modes become mixed. The 
fundamental mode provides first maximum mixedness at t ~ 0.69s (see Fig. 6(a)). Actually, 
at this value of the interaction time, i.e. t ~ 0.69s, the mean-photon numbers for the three 
modes are equal and consequently their purities are equal, i.e. Tr p 2 (t) ~ 0.37 (we have 



20 



checked these facts). When the interaction proceeds, switching energy between waveguides 
continues until tfi = n/2 then the energies in the waveguides are completely interchanged 
(cf. ( |34l) ). At this moment the fundamental (second) mode is in partial pure (maximum 
mixed) state since its mean-photon number is small (large). This behaviour is repeated 
until fit = 7r where the fundamental and second modes provide their initial mean-photon 
numbers. 

Now we draw the attention to find analytically the values of the interaction time for 
which the modes of the system are completely disentangled, i.e. Trp^(t) = 1. This can be 
realized from ([9]) since all calculations in the paper depend on this expression. Therefore 
the condition for disentanglement is that 

(fl Ap A™> ) = f[ (Ap A™> ) . (35) 
i=i j=i 

From the analysis given in the paper this occurs when xt = mis, m is integer. In this 
case the modes are completely disentangled and each of which has a coherent light with 
time-dependent amplitude. On the other hand, the values of interaction times for which the 
system provides its initial light can be obtained by solving simultaneously the following two 
equations 

fit = 2m'ir, xt — m7r > (36) 

where vn! and m are integers. These equations are connected with the requirements that 
dj(t) = oij. Equations fl36l) lead to the fact that the system tends to the initial stage when 
fi = X an d fit = 2m'ir. In this case the linearity of the system compensates the nonlinearity. 

V. CONCLUSION 

In this paper we have discussed the quantum properties for the three-mode codirectional 
nonlinear Kerr coupler, when the modes are initially prepared in coherent light. After ob- 
taining the exact solution of equations of motion we have investigated single- and two-mode 
quadrature squeezing, W function and purity. Various interesting effects have been ob- 
tained. Different forms of squeezing have been achieved by varying the parameters x> <\j an d 
atj. Squeezing can be equally shared between the guided modes and the interacting modes 
can behave as two separated dynamical systems by controlling the intensity of switching be- 
tween waveguides. Quadrature squeezing can exhibit leaf-revival-collapse phenomenon based 
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on the competition between the linearity and nonlinearity in the system. Moreover, different 
forms of cat states have been generated and confirmed in the evolution of W function. We 
have analytically proved that two different types of cat states can be simultaneously realized 
in the evolution of the fundamental and second (third) modes. Also we have discussed the 
mixedness for the single-mode case and deduced the values of the interaction time as well as 
the conditions required for complete disentanglement between the components of the system. 
The nature of the coupler has manifested itself in all studied quantities as complementary 
between the behaviours of the fundamental and second-third modes, i.e. as switching energy 
between waveguides. Moreover, we have showed that the nonclassical effects provided by 
the fundamental mode are richer than those for both the second and third modes as well as 
for the two-mode KNC Finally, the system discussed in this paper is more effective than 
the conventional Kerr coupler and can be used to amplify the nonclassical effects. 
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